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^ Abstract 

In the heavy-quark limit, the decay spectrum for radiative and semileptonic in- 
clusive .B-meson decays is determined by a universal structure function, i.e. the shape 
^ ■ function. We study the constraints on the bilocal heavy-quark operator for the shape 

function from the QCD equations of motion and heavy-quark symmetry, and we obtain 
a new basis of nonlocal operators that reveals relevant soft components in the shape 
function. Those nonlocal operators represent the "kinetic energy distribution" of the b- 
quark inside the 5-meson and the four-parton correlations with additional quarks and 
gluons. A corresponding local operator basis relevant for relating the shape function 
to HQET parameters is also given, and novel effects of the quark-gluon correlations 
are discussed. 
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§1. Introduction 



Inclusive £>-meson decays, such as charmless semileptonic B — > X u £u decays and penguin- 
induced B — > X s 7 decays, are of special interest because they are sensitive probes of elec- 
troweak parameters as well as new physics. Theoretical analysis of the corresponding (dif- 
ferential) decay rates has been developed on the basis of the operator product expansion 
(OPE) as a power series in Aq CD / m b .^ This provides justification for the parton model and 
also the power corrections to it in terms of matrix elements of local operators. 2 ) It is known, 
however, that the usual OPE becomes singular near the kinematic endpoint E/E ms , x ~ 1 
for the lepton (photon) energy spectrum in the decay B — > X u tv [B — > X s j); the endpoint 
region plays an important role in experimental analysis, especially for the precise determi- 
nation of the CKM matrix element \V u b\ from B — > X u £u, avoiding large backgrounds from 
the decays into charmed particles, B — > XJ,u. In the endpoint region, the hadronic decay 
products evolving from the w-quark (s-quark) in B — > X u tv [B — > X s j) have large energy 
but small invariant mass, so that the hadronic final state becomes jet-like with collinear 
interactions of an outgoing light quark. This implies that the short- distance expansion is 
not applicable, receiving the light-cone singularity, and the OPE has to be reorganized with 
the resummation of the most singular terms. 3 )' 4 -* This is accomplished through the light- 
cone expansion, which is similar to the treatment in deep inelastic lepton-nucleon scattering 
(DIS). As the result, at the leading power of A^cd/"^ the shape of the decay spectrum is 
described by an analogue of the leading twist in the DIS, i.e., by the factorization formula, 
where a structure function corresponding to the nonperturbative, long-distance (~ 1/tIqcd) 
contribution is convoluted with the perturbatively calculable function. 3 )~ 5 ) The structure 
function, called the shape function, is process independent and universal in the sense that 
the same shape function determines the decay spectrum of both B — > X u tv and B — > X s ^. 
More rigorous treatments of the factorization with three relevant mass scales, hard (m&), 
hard-collinear (^m b AQ C D), and soft (/Iqcd), have been presented recently in the frame- 
work of soft-collinear effective theory (SCET). 6 )~ 9 ) Those treatments have shown that the 
perturbatively calculable function is composed of two parts, the hard function due to the 
hard corrections and the jet function due to hard-collinear fluctuations associated with the 
light-quark jet, and that the resulting factorization formulae for B — > X u tv and B — > X s ^ 
are valid to all orders in a s . 5 )~ 8 ) 

One-loop corrections 10 ^ as well as the resummation of large (leading and next-to-leading) 
endpoint logarithms 5 )' 7 )~ 9 )' n )~ 13 ) have been calculated for perturbatively calculable functions 
in the factorization formulae for B — > X u tv and B — > X s ^. Also for the shape function, 
one-loop renormalization group evolution has been studied by many authors. 7 )- 9 )> 13 H 5 ) Un- 
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fortunately, however, our understanding of the nonperturbative boundary value in the cor- 
responding evolution equation is still poor. There have been some attempts to constrain 
the boundary value of the shape function by estimating the relevant nonperturbative effects, 
e.g., trying to relate the moments of the shape function with the fundamental parameters 
in the heavy-quark effective theory (HQET), 3 )' 4 )' 16 )' 17 ) and examining the ambiguity ("in- 
frared renormalons" ) of the perturbation series for the shape function. 5 )> 13 ) Guided by such 
constraints, one would construct ansatze for the shape function, and eventually fit them to 
experimental data to fix the remaining uncertainties. 18 -*' 19 - ) 

In attempting to reach this goal, it is desirable to analyze the operator structure of the 
shape function in QCD and clarify the maximal model-independent constraints among the 
matrix elements of the relevant operators. However, a systematic study for this purpose has 
not yet been carried out. In this paper, we discuss a systematic operator analysis for the 
shape function. Our approach is based on the exact identity for the bilocal heavy-quark 
operator, which we derive from the QCD equations of motion and heavy-quark symmetry. 
Combined with the light-cone expansion of the bilocal operator to separate the longitudinal- 
momentum dependence from transverse one, it is possible to derive a differential equation 
that controls the longitudinal-momentum dependence of the shape function. We find that 
the corresponding differential equation involves the "source terms" given by matrix elements 
of a finite set of higher- dimensional nonlocal operators. Those nonlocal operators represent 
the "kinetic energy distribution" of the 6-quark inside the S-meson and the four-parton 
correlations with additional quarks and gluons. Solving this differential equation, we ex- 
press the shape function in terms of those nonlocal operators that explicitly represent the 
nonperturbative effects relevant for the shape function. 

When deriving the differential equation from the exact operator identity, and also when 
solving the differential equation, we need to specify the boundary conditions for the shape 
function. In this paper, we employ the boundary conditions that the (non-negative) moments 
of the shape function as a function of the longitudinal momentum are finite. These conditions 
are equivalent to assuming that the nonlocal operators relevant to the shape function are 
the generating functions of the corresponding local operators via the Taylor expansion, and 
hold trivially when ignoring renormalization effects. *) However, recent results for the one- 
loop renormalization of the shape function in the MS scheme claim that the renormalization 
effects completely modify the mathematical properties as well as the physical interpretation 
of the shape function. In particular, a "radiative tail" is generated for the momentum 
representation of the shape function, so that all non-negative moments become ultra-violet 

*) It is well known that similar conditions arc satisfied by the structure functions in the DIS, and that 
they are preserved by including the radiative corrections. 23 ^ 26 ^ 
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(UV) divergent. 7 -*' 8 )> 13 )' 14 ) Because of this, the differential equation for the shape function 
and its solution discussed in this paper are exact up to the strong renormalization effects 
of order a s . However, it is important to emphasize that actually our results are useful 
and indispensable beyond the lowest order in the perturbative effects, i.e., for obtaining 
the total behavior of the renormalized shape function: As noted in Refs. 7), 8), 13), the 
strong UV behavior in the renormalization effects implies that, for a consistent treatment 
and interpretation, the renormalized shape function should be further "factorized" into the 
"hard components" involving the radiative corrections and the "soft components" involving 
all the nonperturbative effects. In the matching calculation to perform this factorization, 
the first step is to list a basis of operators to describe the relevant soft components. The 
solution given in the present paper provides a basis of nonlocal operators as well as the 
corresponding local operator basis, including the four-parton correlation operators that were 
previously unknown. Furthermore, our solution provides the complete tree-level result for 
the hard components. 



The remainder of this paper is organized as follows. In §2 we demonstrate that the 
equations of motion and heavy-quark symmetry allow us to express the bilocal heavy-quark 
operator for the shape function in terms of a finite set of higher-dimensional nonlocal opera- 
tors, which provides a new basis of nonlocal operators representing relevant soft components 
in the shape function. We also derive a corresponding local operator basis in §3, which is 
relevant when relating the soft components of the shape function to the HQET parameters. 
In contrast to the previous works, 3 -' -5 -'' 13 )' 16 )' 17 ) our local operator basis is composed of an 
increasing number of operators with increasing dimension; in addition to a few first operators 
with low dimension, whose matrix elements give the well-known HQET parameters, 3 )' 4 )' 16 ) 
we find "new" operators with higher dimensions, whose matrix elements provide the "gen- 
eralized" HQET parameters that represent the "Fermi motion" of the 6-quark as well as the 
four-parton correlation effects inside the 5-meson. We also discuss simple estimates of some 
of these generalized HQET parameters and novel effects of the quark-gluon correlations. In 
§4 we derive an explicit formula for the momentum representation of the shape function. Al- 
though this formula is subject to modification due to large radiative corrections, it deserves 
consideration, because it expresses the "Fermi motion" of the 6-quark in an explicit analytic 
form and the four-parton correlation effects in an integral representation. In §5 we present 
our conclusions. 
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§2. Nonlocal HQET operators and equations of motion 



The factorization formulae for the decay spectrum in the inclusive S-meson decays are 
proved with the two-step matching, i.e., matching QCD onto SCET at a scale ~ m&, followed 



at the second step by integrating out the final-state light-quark jet associated with the mass 



Here z M denotes a light-like vector, z 2 = 0, t> M is the 4- velocity of the 5-meson (v 2 = 1), and 
t = v-z. We choose the Lorentz frame of the system as z^ = (0, z~ , 0±) and v M = (v + , v~ , 0j_), 
corresponding to the case in which the final-state jet momentum points in the "— " direction 
on the light-cone. h v (x) denotes the effective 6-quark field, b(x) ~ exp(— im b v ■ x)h v (x), and 
it is subject to the on-shell constraint fh v = h v . 20 ^' 21 ' Here 



is the path-ordered gauge factor along the straight line connecting the points z and 0. 
For brevity, in the following, we do not show the path-ordered gauge factors connecting 
the constituent fields. We employ a mass-independent normalization of the 5-meson state 
\B(v)}, such that (B(v)\B(v')) = v°(27i/mB) 3 S^(v — v'), with mg the I?-meson mass. The 
functions f(t) and f{oj) give the coordinate and (residual) momentum representations of the 
shape function, respectively. Due to the heavy-quark spin symmetry, there appears no other 
independent function by taking matrix elements of the bilocal operators with other Dirac 
matrices inserted. 

Actually, the shape function (j2-l)l depends on the scale /i at which the nonlocal light-cone 
operator is renormalized (a/^^qcd ^ H ^ ^qcd)- It is known that the radiative correc- 
tions induce a Sudakov-type strong scale dependence on (j2-l|) . governed by the corresponding 
cusp anomalous dimension, 5 )> 13 M 5 ) and, corresponding to this effect, fit) in the MS scheme 
becomes singular for short distances: /(0), as well as all derivatives of f(t) at t — 0, diverge. 
As discussed in §1, however, such singular UV behavior can be factorized from the "soft 
components" of the shape function by an additional matching procedure, and we can ignore 
the renormalization in order to disentangle the operator structure relevant for the soft com- 
ponents. Therefore, in the following, we work with "naive" mathematical properties for ()2-l|) 
that are valid at lowest order in the perturbative effects, e.g., that /(0) = 1, and that fit) 



by matching SCET onto HQET at a scale ~ y/m&ylQCD- 5 ^ 9 ^ The shape function is introduced 





(2-2) 
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can be Taylor expanded about t = 0. With these boundary conditions, the physical meaning 
of f(u) is that it is the distribution of the residual momentum k + = uov + of the heavy quark 
inside the S-meson.*^ A similar approach has been employed in Refs. 5) and 13) to analyze a 
few first local operators in the Taylor expansion of h v (z)h v (0) about = 0, considering the 
small t expansion of f(t). Our approach corresponds to an extension of the analysis given in 
Refs. 5) and 13) by treating the relevant nonlocal operator h v (z)h v (0) directly. Although this 
is formally equivalent to treating all local operators in the Taylor expansion of h v (z)h v (0) 
simultaneously, our nonlocal operator approach has the advantage of making the operator 
structure for the soft components of the shape function most transparent. 

We can utilize the nonlocal operator technique, which has been developed for analyzing 
the (higher twist) nucleon structure functions 23 ^ 26 -* (see also Refs. 28) and 29)*- 1 ). The 
nonperturbative dynamics of the 6-quark, surrounded by the light quarks, antiquarks, and 
gluons, reveals itself as the response of the nonlocal operator h v (z)h v (0) to the change of the 
interquark separation and/or total translation, but the total translation is irrelevant for the 
forward matrix elements. The relevant response of the nonlocal operator is described by the 
exact operator identity 

d - - , f 1 - 
; h v (x)h v (0) = h v (x)v ■ D h v (0) + i / duu h v (x)gG fM/ (ux)v fl x u h v (0) , (2-3) 



where x^ is not restricted on the light cone. D^d^ +igA^ and — — igA^ are the 
covariant derivatives, and G^ v = (i/g)[D^, D u ] is the gluon field strength tensor. Taking the 
B- meson matrix element of the relation (J2-3)) . the first term on the RHS vanishes, due to 
the HQET equations of motion h v v ■ D = 0. We take the light-cone limit, x^ — > z^; for the 
calculation of the LHS, we actually need to extend the definition ()2-lj) to the case in which 
the interquark separation is not light-like, where we have 

(B{v)\h v (x)h v {0)\B(v)) = f(vx) +x 2 F(vx) + C(x 4 ) , (24) 

while, for the second term on the RHS, we introduce a three-parton correlation function 
R(t, u) as 

(B(v)\h v (x)gG^(ux)x v h v (0)\B(v)) = l(x^-— v„) \k(t, u) + 0{x 2 )] . (2-5) 

2 \ v ■ x J L J 

Then we obtain the following constraint equation due to the equations of motion: 

+ 2tF(t) = -t I du u R(t, u) . (2-6) 
at 2 J Q 

*-* Actually, the validity of this "probabilistic interpretation" appears to be critical, due to nonperturba- 
tive as well as perturbative effects. (Sec the discussion in §§4 and 5 below.) 

*) A similar extension has been employed for the B-meson light-cone wavefunctions. 27 ^ 
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This equation is formally similar to the corresponding differential equations for the twist-3 
nucleon structure functions, 26 -* but an important difference is the participation of F(t) in 
()2-6j) . which expresses the effect due to the deviation from the light-cone, as in ()2-4|) . For 
the nucleon case, contributions of this type correspond to twist-4 effects and decouple from 
the equations at the twist-3 level, because the twist defined as "dimension minus spin" of 
the relevant operators is a good quantum number. Although the twist counting of the shape 
function (|2-lj) would be twist-3, such a "conventional" twist is no longer a useful concept in 
the HQET. 

It is straightforward to see that the insertion of an arbitrary Dirac matrix r into ()2-3j) - 
()2-4|) does not lead to a new equation for f(t). This is unlike the case for the i?-meson 
wavefunctions, 27 ) and it is due to so strong constraints on the forward matrix elements 
imposed by the heavy-quark spin symmetry. 

To proceed, we analyze the explicit operator structure of F(t) in (|2-fi|) . which can be 
extracted from the next-to-leading term in the light-cone expansion of the nonlocal operator 
h v (x)h v (0) [see (|2-4|) ]. For this purpose, we utilize an elegant method to construct the light- 
cone expansion: 23 -*' 24 )< 29 ) The leading term (It) in the light-cone expansion, corresponding to 
the leading twist operator in the DIS, obeys the equation 

[h v (x)h v (0)] lt = , (2-7) 



dx^dx^ 



which ensures that all local operators arising in the Taylor expansion are traceless. A formal 
solution is 23 ^ 

[h v (x)h v (0)] h = h v (x)h v (0) 

+£ ^y. jf ^^iP- (a^y " M-)M0) . (2-8) 

To order x 2 , we obtain 



— — P u1l cP" - 

h v {x)h v {0) = [h v (x)h v (0)] lt + -r / — T—r—K{ux)h v {0) + 0{x A ) . (2-9) 
L J " 4 Jo « ox^ax^ 

Actually, the first term on the RHS of ()2-4|) contains the 0(x 2 ) term as well as the leading 
term. To separate contributions of different powers in x 2 , it is convenient to exploit the light- 
cone expansion of exp(iuv ■ x) entering into the definition of the momentum representation 
of />-x), 24 )> 29 ) 

■. = [ e *«-*] _ / duue im + 0(x 4 ) , (2-10) 



e 
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where [e to ^' x ] lt is defined by a straightforward generalization of the procedure [. . .] lt of (J2-8)) 
to an arbitrary function of x. Substituting (|2-9|) - ()2-10|) into ()2-4|) and comparing both sides 
of the resulting equation, the leading terms reproduce the definition ()2-l|) . and we get, from 
the next-to- leading 0(x 2 ) terms, 



F(t) = \f duu<P(ut) + ^ J duouo 2 f{uo) J duue iuujt , (2-11) 



with 

$(t)= (B(v)\^^h v (x)h v (0)\B(v)) 



(2-12) 



Here, the second term on the RHS of (|2T1|) is the analogue of Nachtmann's correction in the 
DIS. An important point is that this term depends on f(u), and it produces an additional 
term involving f(t) in (|2-6j) . Substituting (|2-11|) into (|2-6|) . we obtain 

t-^- + f(t) - 1 = t 2 duu (iR{t, u) - S{utfj . (2- 13) 

Here the RHS involves $ of ()2-12|) . which is concerned with the effects of the transverse 
as well as longitudinal motion of the 6-quark inside the 5-mesons. To reveal its physical 
content further, the following exact identity is useful:*) 

- — — h v (x)h v (0) = h v (x)(Du,) 2 h v (0) + 2i duu- — h v (x) gG^ux) x u h v (Q) 
dx^dx^ p J dx^ 

— if duu 2 h v (x) [D p , gG ^(ux)} x v h v (0) 



o 



+ 2/ duu dssh v (x) gG ilv {ux)x v gG^ p (sx) x p h v (0) . (2-14) 



This can be derived straightforwardly by a method similar to (|2-3j) . Substituting ()2-14|) into 
(|2-12j) . the replacement D^D±^=D P —v^v- D can be made in the first term on the RHS 
by using the equations of motion for the effective heavy-quark field. Similarly, the equations 
of motion for the gluons imply [D^ 1 , G^ u ] = —gt a (v u h v t a h v + qt a ^ v q) , with the summation 
(S) over the heavy (h v ) and light (q) flavors. Therefore, ()2-12|) can be expressed in terms 
of a set of many-body nonlocal operators. This pattern is similar to that for the twist-4 
effects in the DIS, 23 )' 24 - 1 but one important difference is that the covariant derivative for the 
transverse direction, Dj_ p , acting on the quark fields cannot be completely eliminated in the 
present case; participation of such a transverse derivative in a complete set of the higher- 
dimensional operators is typical of the HQET. 20 )' 21 ) In fact, the corresponding operator in 
*> A similar identity for the case of the light quarks is discussed in Ref. 29). 
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(j2-14j) plays an important role in describing the Fermi motion effects in the 5-mesons, as we 
show below. 

The matrix element of the derivative of the three-body operator in (j2-14|) can be calcu- 
lated using (J2-5)) . An important observation is that this contribution cancels the first term 
iR(t,u) on the RHS of (j2-13j) . so that the three-parton correlation R decouples from our 
differential equation. For the remaining terms, we define 

(B(v)\h v (z)(D ±fl ) 2 h v (0)\B(v)) = (B(v)\h v (z)(D^) 2 h v (Q)\B(v)) =K(t) , (2-15) 
(B{v)\h v {z) {D»,gG^(uz)]z»h v (0)\B(v)) 
= -J2(B(v)\hv(z)g 2 t a q(uz)t a ^q(uz)h v (0)\B(v)} =tW(t,u) , (2-16) 

(B{v)\h v {z)gG^{uz)z v gGW{sz)z^)\B{v)) = t 2 Y{t,u,s) . (2-17) 

Here, K(t) can be interpreted as the "kinetic energy distribution" of the 6-quark inside the 
S-meson, while W(t, u) and Y (t, u, s) are the four-parton correlations with additional quarks 
and gluons, respectively. We denote the RHS of (|2-13j) as J{t) and substitute (|2-15|) - (j2-17l) 
into it. We then obtain 

J{t)= [ dt' l-t'kit') +it' 2 [ duu 2 W{t',u)-2t' 3 [ duu [ dss Y(t',u,s) \ . (2-18) 
Jo I Jo Jo Jo J 

Regarding J(t) as the "source" term, (j2-13J) is immediately integrated to give, with the 

condition f(0) = 1, 

f{t) = l + \ I drJ(r) . (2-19) 



tjo 

Thus f{t) is completely re-expressed in terms of a set of nonperturbative functions, K(t), W(t, u) 
and Y(t,u, s), corresponding to the higher- dimensional nonlocal operators. 

§3. Local operator basis and HQET parameters 

In this section we derive a basis of local composite operators to represent the soft com- 
ponents of the shape function, and we discuss the HQET parameters as matrix elements of 
these operators. It is straightforward to expand the solution ()2-19|) in a power series in t as 

L — ' n! 

n=l 

where the n-th term gives the nonperturbative power corrections of order (/Iqcd^)™- We 
obtain 



a ~ n ~ 1 r 1 



n + 1 n(n + 1) 



J2 {(* + !)(* + 2) 



k=0 
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+ (n — k — l)(n — k - 2)}W k , n -3-k - 



l) n n 2 - 1. 



W n-3 n-3 

4 2 ' 2 



ra(n + 1) 

i + (-r 



n— 4 

E 

k=0 

k / k 



^2 {(l + l)(n + 1 — k — 1) + (k — l + l)(n — l — 1)} y t 

1=0 

k 



n—i—k 
Lk-l 



+ 1 



n - - - 1 K* 

Z / 2 ' 2 



4-fe 



(3-2) 



for n = 1,2,3,.... Here /C n ,W/ > fc and 3^™ denote matrix elements of the local operators, 
which are generated from the Taylor expansion of (|2-15|) - (|2-17|) about = 0: 



= j^(B(v)\h v (tD + ) n (D ± ,) 2 h v \B(v)) , 
m,k = —, + * +1 ^(B(v)\h(iD + )Wt a &7MiD + )%\B(v)) 



^i,k 



■(B(v)\h v (iD + ) k gG, + (iD + rgG^(iD + yh v \B(v)) . 



(3-3) 
(3-4) 

(3-5) 



( v +^k+l+m+2 

Note that Wi t k and 3^™ are symmetric under the interchange I <-> k, which follows from the 
behavior of the corresponding matrix elements ()3-4|) and ()3-5|) under the parity transforma- 
tion combined with the time- reversal transformation. Equation ()3-2j) gives, for n = 1 and 2, 



Ai = 0, 



•Al = — /Cn = - 

3 3 



(3-6) 



with the fundamental HQET parameter Ai related to the average kinetic energy of the heavy 
quark inside the S-mesons, 20 -'' 21 '' 



and, for n = 3, 



Ai = (B(v)\h v ^D ±fl ) 2 h v \B{v)) 



A 3 = -i>Vo,o = \Y,(B{v)\h v gh a qt a i>qh v \B{v)) , 
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(3-7) 



(3-8) 



where the four-quark operator is related to the quark-gluon Darwin operator by the equations 
of motion [see also (|2-16|) ]. Equation (|3-6|) coincides with that obtained in Refs. 5) and 13). 
Also, ()3-8|) as well as ()3-6|) coincides with the corresponding results discussed in the earlier 
works. 3 -'' 4 )' 16 ) We note that those works analyzed the constraints resulting from the equations 
of motion on the matrix elements of the corresponding local operator, 

1 



A 



(V+) 



-{B(v)\h v (iD + ) n h v \B{v)) 



(3-9) 



which follows from ()2-l|) . 
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Furthermore, ()3-2j) gives new results for n > 4, 

A,= \K 2 - 2 ~W , x -^y%, (3-10) 

As = ~ TE m ' 2 ~ ~ l^ ~ I^ 1 ' (3,11) 

and so forth. Therefore, the operators in (|3-3jl - ()3-5jl . whose matrix elements give 

n — 3 



/C n _2, W k ,n~3-k [k = 0, 1, 2 



2 



y^if* ^ = 0,1,2,- -.,71-4; / = 0,1,2,---, |) 



(3-12) 



form a basis of local operators of dimension n + 3 and contribute to ^4 n . This implies that 
for n ^> 1 there are ~ n 2 /4 independent local operators participating in A n . 

Here it is worth studying the details of the new nonperturbative matrix elements con- 
tributing to Aa, because these matrix elements provide a generalization of the HQET pa- 
rameters. K.2 and Wo,i, given by ()3-3|) and (|3-4j) . can be related to the "trace part" of the 
corresponding operators as 

Kz = \{B{v)\h v ({D^ff h v \B(v)) , (3-13) 

W 0)1 = l - J2( B ( v )\^9 2 t a qt a ^qD^h v \B(v)) . (3-14) 
g 

Thus JC2 > 0, and a comparison of (j3-13|) with ()3-7|) leads to the rough estimate K.2 ~ A 2 . 
The matrix element ()3-14j) involving the four-quark operator can be evaluated in the vacuum 
saturation ("factorization") approximation*) as 

N 2 -1 

W ,i « -7T ^ 2 a 8 f 2 B m B A , (3-15) 

with the usual decay constant defined as (0\q , -f IJ, 'y 5 h v \B(v )) = ifB\/m B /2 v^. Here j4 
represents the asymptotic value of the mass difference mg — mj between the I?-meson and 
the 6-quark, 20 -'' 21 '' and it can be identified with the effective mass of the light degrees of 
freedom in the 5-meson. Using phenomenological values for the parameters on the RHS of 
(j3-15|) leads to Wo,i ~ — A 2 . Contrastingly, 3^o,o °f P'^P cannot be estimated in a simple 
manner, but in the 5-meson rest frame it is expressed in terms of the chromoelectric and 
chromomagnetic fields as 

y°o,o = -\{B{v = 0)\h v g 2 (E 2 + B 2 ) h v \B(v = 0)> , (3-16) 



*) This approximation gives the estimate W ,o ~ - l)/2N^]a s f^m B of (|3-8|) 



4), 16) 
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so that 3^o,o < 0- Here we note that a comparison between (j3-6j) and (J3-9)) with n = 4 
suggests the rough estimate .4.4 ~ Af, and combined with the relation (|3-1U|) .*) we would 
expect 3^o,o ~ ~~ 111 principle, A5 of (|3-11|) could be analyzed similarly, but it involves the 
five new nonperturbative matrix elements and is much more complicated. 

When AqcDt is not so small compared to 1, which is the case in the endpoint region 
for the lepton (photon) energy spectrum in B — > X u lv (B — > X a j), the higher-order power 
corrections {AqcDt) 71 are important**-* to determine the behavior of /(£). Among the param- 
eters ()3-3|) - ([3-5|) involved in A n , the experimental value is known for /C (= — Ai) of (|3-6|) 
only , 8 )' 31 ) Hopefully, further information regarding these generalized HQET parameters will 
be obtained from nonperturbative calculations with, e.g., lattice QCD***) and QCD sum 
rules. 

§4. Fermi motion and four-parton correlation effects 

In this section we derive an explicit formula for the momentum representation of the 
shape function, which is exact up to perturbative 0(a s ) corrections. Although this formula 
suffers from (large) perturbative corrections, it is useful for obtaining some insight into 
the relevant nonperturbative effects, the Fermi motion of the 6-quark, and the four-parton 
correlation effects, without recourse to the t-expansion as in (|3-1|) . For this purpose, we 
derive the Fourier transformation of ()2-19|) into the residual momentum representation [see 
(J2HJ)]- We obtain f{ui) = 5(u) + [contribution from J]. Clearly, the first term, S(u), reflects 
the lack of heavy-quark motion inside the 5-mesons in the "non-interacting limit", and 
the above decomposition of f(ui) into the "free" term S(u) and the "interaction-dependent" 
contribution from J would be unsuitable to treat the nonperturbative part in the shape 
function. In order to avoid such an "artificial" singularity, we rearrange the terms in (|2-19j) . 
Actually, it is more convenient to make the corresponding rearrangement in our differential 
equation (|2-13J) by moving a certain term from the "source" on the RHS to the LHS, and 
then solve the resulting equation with the new source. We decompose K (t) in (j2-18j) as [see 

Kit) = a 2 fit) + 6K(t) , (4-1) 

*) The identity between matrix elements of the local operators, which is implied by l|3-9|l with n = 4 
and (|3 • 10|) . would be subject to the perturbative 0(a s ) corrections, but these corrections do not affect this 
simple estimate. 

**J For the regions with m& ^> 1/t, these power corrections are enhanced compared with the usual sublead- 
ing corrections, which are suppressed by powers of Aqqb / Tib- 32 -* Physically, this corresponds to the situation 
in which the decay spectra in the endpoint region is smeared over a range A with m& 3> A > ^Iqcd- 7 ^' 9 ^ 
***) p or recent progress in treating matrix elements of B-mesons in lattice QCD, see, e.g., Ref. 33). 
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and move the first term with the constant a 2 to the LHS of ()2-13j) . The residual term 8K(t) 
is treated as a source term. The explicit form of the residual term reads 

5K(t) = (B(v)\h v (z) ((D^) 2 - a 2 ) ^(0)|S(t;)) . (4-2) 

The rearrangement using (|4-ljl is motivated by the observation that the first term is special 
among the terms on the RHS of (J2-18|) : Among the nonlocal operators of (|2-15J1 - (|2-17J1 
contributing to (|2-18|) . the operator for K(t) has the lowest mass dimension, and it is expected 
to play a dominant role. In fact, only this contribution is of a bilocal nature as (j2-15|) . as 
in the case of the shape-function operator ([2-ljh while the other contributions come from 
genuine many-body operators. Noting the properties of the higher twist nucleon structure 
functions in the DIS, 30 ^ we conjecture that multi-parton correlation effects with increasing 
numbers of partons are less important. We note that the operator for K(t) is a nonlocal 
version of the kinetic energy operator of (|3-7j) . and it should be relevant to the Fermi motion 
of the 6-quark inside the 5-meson; in (J4-lj) we have extracted the part proportional to f(t) 
from K(t) on the basis of the fact that K (t) and f(t) have similar bilocal natures, and thus 
a would eventually correspond to the "average measure" of the transverse quark motion, 
as a 2 ~ — Ai [note that -^(0) = — Ai]. However, for the time being, we can leave a as an 
arbitrary parameter. 

In momentum space, our differential equation now reads 

{u 2 - a 2 ) = _/ M . (4-3) 

Here, the source I(co), defined as dJ(t)/dt _ =i J duje lu)t I{uj), is given explicitly as [see 

GEEK] 



d r a r°° i ( d 2 d 2 \ 

m = --sk M + / da p / „ _ (_ + _) 



W(a,P) 



+ 2 daP df3 d£ 



a — (3 



1 f d_ d_ 



a — £ \da <9£ 



d_ d_Y 1 / d d_ 



1 f d d 

+ 7 Zw"7i ^ 7^ + 7^: 



(a- W -0 \dP 1 dU1 Y ^M-T 2 l Y{a -^ M ] ' (44) 

where "P" denotes the principal value, and we have introduced the momentum representation 
for the relevant matrix elements ()4-2|) . (|2Tfij) . and (|2-17j) as 

6K(t) = [ du e iwt 5K{u) , (4-5) 
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W(t,u) = J duj (ku'e^-^'-^Wiu,^) , (4-6) 
Y(t,u,s) = [ dujduj'd^ e^ t -W-® u *- i{ £- u) ' t Y(u,u',£) . (4-7) 



We note that the parity transformation combined with the time-reversal transformation 
implies the symmetry properties 

W(u, J) = W(uj', u), Y(uj, uj', = Y(uj', uj, f) , (4-8) 

and hermiticity guarantees that W(uj,uj') and Y(u,uj',£) are real functions. For both func- 
tions, the variable uj has the physical meaning of the residual light-cone momentum carried 
by the 6-quark. The quantity uj' — uj for W(uj,uj') represents the light-cone momentum car- 
ried by the gluon or the light gg-pair, while £ — uj and uj' — £ for Y(uj,uj',C,) represent the 
light-cone momenta carried by the gluons. By inserting a complete set of states between 
the constituent fields, we also get the support property: W(uj, uj'), as well as Y(u,w',£), 
vanishes unless uj < A and uj' < A. Similarly, we deduce that 8K{u) is a real function that 
vanishes unless uj < A. Using these symmetry and support properties of SK, W and Y, it is 
straightforward to see that I{uS) of ()4-4|) vanishes unless uj < A. 

We solve ()4-3|) with the boundary conditions f{uj) = for uj — > ±oo and the normalization 
condition f dujfiyj) = /(0) = 1 in the form 

f(uj) = f^(uj) + fW(uj) , (4-9) 

where f(°\uj) is the solution with the source I(uj) set to zero, while f^\uj) denotes the piece 
induced by I(uj). It is straightforward to obtain the analytic solution 

/<W( W ) = ±.{0(u> + a)-0{u>-a)} , (4-10) 
la 

which satisfies J dujf^\uj) = 1. Also, the solution for f^(uj) reads 

/«(„) = ± ( r~' ] *j ^ - r + " ^^\, (4.11) 



2cr yj u uj' -a J,, uj' + a 

which satisfies J du (uj) = 0. The solution ()4-9|) with (|4-10|l and (|4T1J) represents an 
explicit formula for the shape function, which satisfies the HQET equations of motion exactly. 
We note that, if we take the limit a -> 0, we get f i0) (uj) -> 5(uj), so that (gJD|) and (l4~TTll 
would reduce to the "singular" decomposition of f(uj) without the rearrangement using (|4-lj) . 
Therefore, ()4T0|) involves a nonperturbative effect represented by the parameter a, which 
smears the singular behavior of 5(u). If we choose, e.g., a = y/—Xi, the corresponding effect 
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would be interpreted as the Fermi motion effect of the 6-quark inside the 5-meson. This 
choice can be useful in practice, because, in the coordinate representation, we have 

and the expansion on the RHS reproduces the first three terms of (j3-l|) exactly [see (|3-6jl . 
and note that /«(t) = it 3 W , /36 + 0{t% In fact, "flat distribution" (14"TUD with 

cr = \/— \\ exhibits behavior consistent with the Gaussian distribution around uo = with 



the width y— Ai, which was proposed in Refs. 5) and 13) as a nonperturbative ansatz valid 
in the vicinity of u = 0. As an alternative to the choice a = y/—Xi, one may employ another 
"optimized" choice for a, which makes 5K(t) of (|4-2|) "small", so that 5K(t) could be treated 
perturbation. 

It is well known that ()2T|) should obey f(uj) = for uj > A 3 )' 4 )' 16 ) Using the above- 
mentioned support property for I(u), it is straightforward to show that our solution ()4-9|) 
indeed vanishes unless to < A*^ reproducing the correct support property. 

Our solution reveals that the shape function is subject to nonperturbative effects due 
to the multi-parton correlation with additional quarks and gluons, which are represented 
by W(cu,uj') and Y(cj,u',^), as in (|4-11|) . (|4-4|) . This indicates that f(u>) is actually a 
complicated multi-particle object beyond a "simple" momentum distribution function. This 
is in contrast to the leading-twist nucleon structure functions in the DIS, but it is reminiscent 
of the character of the higher-twist ones. 23 ^ 25 ^ 26 - 1 Our results, (|4-llj) and (|4-4jl . are exact 
to lowest order in the perturbative effects, but the multi-particle character of f(ou) should 
persist when perturbative corrections are included. 

§5. Conclusion 

The shape function of the 5-meson is an important ingredient in the factorization for- 
mula of the differential rates of inclusive 5-meson decays, such as B — > X s 7 and B — ► X U W. 
We studied the "soft components" of the shape function separately from the "hard (pertur- 
bative) components" relevant to its scale dependence, and we studied the model-independent 
constraints on the soft components in the heavy quark limit. From the equations of motion 
and heavy-quark symmetry constraints, we derived a differential equation that relates the 
longitudinal-momentum dependence of the shape function to matrix elements of the novel 
nonlocal operators. Solving this equation, we obtained a new basis of nonlocal operators, 
which describes the relevant soft components of the shape function in terms of the kinetic 



*' This is actually satisfied even in the case a > A, as well as for a < A. 
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energy distribution and the four-parton correlations. We also derived the relations among 
the matrix elements of local operators with arbitrary numbers of covariant derivatives. These 
relations not only reproduce all the known results for a few covariant derivatives but also 
include new relations for more covariant derivatives. These relations are exact only when 
we ignore the perturbative effects, but they are useful beyond the lowest order in the per- 
turbative effects for the purpose of finding a basis of independent local operators with the 
same dimension. Our relations indicate that the local operator basis is composed of an 
increasing number of operators with increasing dimensions whose matrix elements give the 
generalized HQET parameters that represent the Fermi motion of the 6-quark as well as the 
four-parton correlation effects inside the 5-meson. Our differential equation also yields the 
momentum representation f(u) as a sum of the part f(°>(uj) involving the Fermi motion and 
the additional part f^(u) induced by more sophisticated correlation effects. The behavior 
of the "leading" part, f(°\u), appears to be consistent with the existing nonperturbative 
ansatz 5 )' 13 ) for the shape function. Our momentum-representation formula is again exact 
up to perturbative corrections, but it provides insight into model building of the relevant 
nonperturbative effects. 

Our results reveal that, due to nonperturbative effects, f(u) is a much more complicated 
object than the simple momentum distribution of the 6-quark inside the 5-meson: Nonper- 
turbative effects induce the mixing of additional dynamical quarks and gluons, so that the 
shape function contains the soft components that represent the multi-particle correlation ef- 
fects. In connection to this, it is interesting to note that a probabilistic interpretation of f(u>) 
as the momentum distribution has been questioned from a different point of view, noting 
that J duif{uj) is negatively divergent, due to the perturbative effects in (j2-lj) . as mentioned 
in §1. 7 )< 8 ) 

The "full" shape function (|2Tj) entering into the factorization formula for the inclusive 
decay rates is obtained by combining the soft components with the hard components through 
the matching procedure. We emphasize that (|3-12|) with (|3-3|) - (j3-5|) provides a local operator 
basis that is necessary for parameterizing the soft components in the matching calculation. 
Furthermore, ()3-l|) with (|3-2|) completes the matching at the tree level for operators of 
arbitrary dimension, and the corresponding leading-order matching coefficients can be easily 
read off from this formula. One can proceed to the one-loop matching calculations using 
our operator basis, and the loop effects will produce the mixing coefficients at the next-to- 
leading order between the operators with the same dimension. Also, the nonlocal operator 
basis provided by (EHHD - (l2~T7J) may be useful, and (l2~T?H) (flUD with (gTXJ and (PETTI) ) 
completes the corresponding leading-order matching in the coordinate (momentum) space. 
Using the results of the matching calculations, a detailed study of the interplay between 
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the strong perturbative effects in the hard components, which generate a Sudakov-type 
scale dependence, and the nonperturbative effects in the soft components, which have been 
unraveled in this paper, should clarify the behavior of the full shape function and its roles 
in decay rates. This will be studied in a separate publication. 
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